Stability and phase coherence of trapped ID Bose gases 
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We discuss stability and phase coherence of ID trapped Bose gases and find that inelastic decay 
processes, such as three-body recombination, are suppressed in the strongly interacting (Tonks- 
Girardeau) and intermediate regimes. This is promising for achieving these regimes with a large 
number of particles. "Fermionization" of the system reduces the phase coherence length, and at 
T = the gas is fully phase coherent only deeply in the weakly interacting (Gross-Pitaevskii) regime. 
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Recent success in creating quantum degenerate ID 
trapped atomic gases 0, 0, U has stimulated an interest 
in correlation properties of these systems. The ID regime 
is reached by tightly confining the radial motion of atoms 
in a cylindrical trap to zero point oscillations. The differ- 
ence between such a kinematically fD gas and a purely 
fD gas is related only to the value of the interparticle 
interaction, which depends on the radial confinement. 

The ID Bose gas with repulsive short-range interac- 
tions characterized by the coupling constant g > ex- 
hibits remarkable properties. Counterintuitively, it be- 
comes more non-ideal with decreasing ID density n 0,0]- 
The equation of state and correlation functions depend 
crucially on the parameter 



7 = mg/h 2 n, 



(1) 



where m is the atom mass. For comparatively large 
n, the parameter 7 < 1 and ^7 represents the ratio 
of the mean interparticle separation 1/n to the corre- 
lation length l c = h/y/mng. In this case the gas is in 
the weakly interacting or Gross-Pitaevskii (GP) regime, 
and the amplitude of the boson field obeys the familiar 
Gross-Pitaevskii equation. For sufficiently low densities 
(or large g), one has 7 » 1 and the strongly interact- 
ing or Tonks- Girardeau (TG) regime is reached. In this 
regime the ID gas acquires fermionic properties in the 
sense that the wave function strongly decreases as par- 
ticles approach each other (see 0, HI and cf. [fj ) . For 
7 = 00 any correlation function of the density can be 
calculated straightforwardly Q by using the exact map- 
ing onto the system of free fermions found by Girardeau 



A uniform ID system of bosons interacting via a delta- 
functional potential is integrable by using the Bethe 
ansatz and has been a subject of extensive theoretical 
studies. Lieb and Liniger [3 have calculated the ground 
state energy and excitation spectrum for any value of the 
parameter 7. Yang and Yang || have studied thermo- 
dynamic properties of this system at finite temperatures 
and found no phase transition for T > 0. The absence 
of long-range order and true Bosc-Einstein condensate in 



finite-temperature ID Bose gases follows from the Bo- 
goliubov fc~ 2 theorem as has been expounded in 0. The 
long-range order is destroyed by long-wave fluctuations 
of the phase leading to an exponential decay of the one- 
particle density matrix at large distances 0, [n|. A 
similar picture is found at T = 0h where the density 
matrix undergoes a power-law decay 0, 0, Q] . In the 
past few decades, a general approach has been developed 
for exact calculations of one- and two-particle correlation 
functions at an arbitrary 7 (see |l5l Il6j for review) . 

Realization of ID trapped Bose gases raises the ques- 
tion of their phase coherence and stability. The phase 
coherence properties are strongly influenced by the trap- 
ping potential, which introduces a finite size of the system 
and provides a low-momentum cut-off of the phase fluc- 
tuations. In the GP regime at sufficiently low T, the 
phase fluctuations are suppressed and the equilibrium 
state is a true condensate. At higher temperatures, it 
is transformed into a phase- fluctuating condensate |17j . 
The dynamics, interference effects, and excitations of ID 
tra ppe d Bose gase s are currently a subject of active stud- 
ies [13, UJi |2(j, |21j . Of particular interest is the change in 
phase coherence properties while going from the GP to 
the TG regime, where the phase coherence is completely 
lost. 

The strong transverse confinement required for the ID 
regime can lead to high 3D densities of a trapped gas. At 
a large number of particles, the 3D density can exceed 
10 15 cm -3 , and one expects a fast decay due to three- 
body recombination. It is then crucial to understand 
how the correlation properties of the gas influence the 
decay rate. 

In this Letter, we discuss stability of ID Bose gases and 
calculate local density correlators, as those are responsi- 
ble for inelastic decay processes 22]. We find that the 
decay rates are suppressed in the TG and intermediate 
regimes, which is promising for achieving these regimes 
with a large number of particles. We then analyze phase 
coherence of a trapped ID Bose gas and show that vac- 
uum fluctuations of the phase make the zero-temperature 
coherence length smaller than the Thomas-Fermi size of 



2 



the sample, unless the gas is deeply in the GP regime. 

The ID regime in a trapped gas is realized if the ampli- 
tude of transverse zero point oscillations Iq — \Jfij mtoo 
is much smaller than the longitudinal correlation length 
l c = fi/y/mjl, where ujq is the frequency of the transverse 
confinement and the chemical potential of the ID system 
is /i <C hujQ. One then has a ID system of bosons inter- 
acting with each other via a short-range potential char- 
acterized by an effective coupling constant g > 0. This 
constant is expressed through the 3D scattering length a 
[23j | , assuming that Iq greatly exceeds the radius of inter- 
action between atoms. For a positive a <C Iq we have 



g = 2h 2 a/mlQ, 



(2) 



and a characteristic distance T% 2 /mg related to the inter- 
action between particles in the described ID problem is 
~ I 2 , /a ;§> Iq. In the GP regime, the chemical potential 
fi w gn, and the condition l c ^> Iq leads to the inequal- 
ity na « 1. In the TG regime the correlation length 
l c ~ 1/n, and one should have jiIq <C 1. We thus see that 
at any value of 7 it is sufficient to satisfy the inequali- 
ties a <C Iq <C 1/n. Then the ID regime is reached and 
the system can be analyzed on the basis of the ID Lieb- 
Liniger model, assuming a delta-functional interatomic 
potential with the coupling constant g given by Eq. J5J. 

The rate of three-body recombination is proportional 
to the local three-particle correlation function 53 = 
{&(x)&(x)&(x)V(x)V(x)V(x)) H2, where *(x) is the 
field operator of atoms and the symbol (. . .) denotes the 
expectation value. Similarly, the rates of two-body in- 
elastic processes involve the correlation function g 2 = 
(^>' (x)^ (x)ty (x)ty (x)) . Assuming that local correlation 
properties are insensitive to the geometry of the system 
we consider a uniform ID gas of N bosons on a ring of 
circumference L. The Hamiltonian of the system reads: 

H = J dx [(h 2 /2m)d x ^d x ^ + (g/2)¥^W] . (3) 

For findingoo at T = 0, we use the Hellmann-Feynman 
theorem |24ll25| . This is similar to the calculation of the 
mean interaction energy in ref. ||. Namely, one shows 
that the expectation value of the four-operator term in 
the Hamiltonian J2J is proportional to the derivative of 
the ground state energy with respect to the coupling con- 
stant: dE Q /dg = ($ \dH/dg\$ Q ) = g 2 L/2. The first 
identity follows from the normalization of the ground 
state wave function <E>o, and the second one is obtained 
straightforwardly from the Hamiltonian (0) . The ground 
state energy can be written as E — Ne(~/)h 2 n 2 /2m, 
where the quantity e{^) is a solution of the Lieb-Liniger 
equations [al and is calculated numerically for any value 
of 7 0, l2fij . For the two-particle local correlation func- 
tion we then obtain 52(7) = n 2 de(p[) / d^ . The function 
92{l)/n 2 is shown in Fig. ^ For small values of 7, we 
obtain numerically the result which coincides with that 
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FIG. 1: Local correlation function (72 versus 7. 



following from the Bogoliubov approach: 

g 2 ( 1 )/n 2 = l-2^/n, 7 « 1. (4) 

In the limit of large 7, we have e(j) — (7r 2 /3) (1 — 4/7), 
and the two-particle correlation function is given by 



g 2 (l)/n 2 = 4^ 2 /3 7 2 , 7 » 1. 



(5) 



The results in Fig. ^ and Eq. (J3J) clearly show that two- 
particle correlations and, hence, the rates of pair inelastic 
processes, are suppressed for 7^1. This provides a pos- 
sibility for identifying the TG and intermediate regimes 
of a trapped ID Bose gas through the measurement of 
photoassociation in pair interatomic collisions. 

The three-particle local density correlator g^ cannot 
be obtained from the Hellmann-Feynman theorem. In 
the weak coupling GP regime (7 « 1) one can use the 
Bogoliubov approach, which immediately yields 



93(l)/n 3 ~ 1 - 6V7/7T, 7<1. 



(6) 



For the TG regime (7 ^> 1), we have developed a 
method for calculating the leading behavior of local den- 
sity correlators. Details will be given elsewhere, and here 
we present a compact derivation of 53 at T = 0. In first 
quantization the expression for this function reads 



33(7) 



dX 



^ ) (0,0,0,x i ,...,x N ) ,(7) 



3!(JV-3)! 

where dX — dx<± . . . dx n , and ^q 7 -* is the ground state 
wave function given in the domain < Xi < . . . < xm < 
L by the Bethe ansatz solution: 

(x 1 ,x 2 ,---,x N ) oc ^a(P)cxp ji^fcp^j j , (8) 
p 

where P is a permutation of ./V numbers, quasimomenta 
kj are solutions of the Bethe ansatz equations, and 



a(P) 



n 

3<l 



i^n - k Pj + k Pl 



For 7 ^> 1, we extract the leading contribution to $[ ) 7 ' 1 at 
three coinciding points by symmetrizing the amplitudes 
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a(P) over the first three elements of the permutation P: 



(1771)" 



(9) 



3<l 



where P p = P Pl , P P2 , P P3 , P 4 , . . . P Nl and j, / = 1,2,3. 
The sign of the permutation P is ep, and p runs over 
six permutations of 1, 2, 3. For large 7, the difference 
of quasi-momenta kj from their values at 7 = 00 is of 
order I/7 and can be neglected. Then, from Eqs. © 
and ©, we conclude that to this level of accuracy the 
ground state wave function at three coinciding points is 
given by derivatives of the wave function of free fermions 



^°°\xi,X2,x 3 ,Xi, . . .) at x\ 



X2 



x 3 = 0: 



$ ( 7) (0,0,0,x 4 ,...) 



( 7 n) ; 



■[nw* 



d Xl ) ^-(lO) 



3<l 



Substituting Eq. (fTUfl into Eq. J7J) we express the lo- 
cal correlator g 3 through derivatives of the three-body 
correlation function of free fermions. Using Wick's 
theorem, the latter is given by a sum of products of 
one-particle fermionic Green's functions G(x — y) = 
/_jt dk e lk ( x ~ v ^ /2tt, where kp = im is the Fermi 
wavevector. The calculation from Eq. (JJJ) is then straight- 
forward and we obtain (7 3> 1) 



53(7) 
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where G and its derivatives are taken at x — y = 0. 

In fact, in our problem the correlation functions gi 
and £3 are slightly nonlocal. They are related to inter- 
particle distances r ~ Z <C ? C ) since at smaller r the 
relative motion of particles is three-dimensional and the 
local correlators do not change. For large 7, at distances 
~ Iq an extra (coordinate-dependent) contribution to 92 
is of the order of (nlo) 2 and to g 3 of the order of (nlo) 6 
(see, e.g. 7J and references therein). Under the condi- 
tion Iq 3> a, these contributions can be neglected as they 
are much smaller than the results of Eqs. (JSJ) and (|llf> . 
respectively. 

Our method is readily generalized for the case of finite 
temperature by considering the temperature-dependent 
Green's functions of free fermions. For T<^we obtain a 
small correction ~ (T//i) 2 to the zero-temperature result. 
The same conclusion holds for the GP regime (7 <C 1). 

Thus, from Eq. (jl 1|> we conclude that the three-body 
decay of ID trapped Bose gases is strongly suppressed in 
the TG regime. Moreover, Eq. © shows that even in the 
GP regime with 7 w 10~ 2 , one has a 20% reduction of 
the three-body rate. Thus, one also expects a significant 
reduction of the three-body decay in the intermediate 
regime. 

For Iq 3> a, the recombination process in ID trapped 
gases occurs at interparticle distances much smaller than 



Iq- Therefore, the equation for the recombination rate 
is the same as in 3D cylindrical Bose-Einstein conden- 
sates with the Gaussian radial density profile. There is 
only an extra reduction by a factor of g 3 /n 3 . A char- 
acteristic decay time r is then given by the relation 
t _1 = a 3 Dn\ D (g 3 / in 3 ) , where a 3 n is the recombination 
rate constant for a 3D condensate, and n 3 n — n /(^o) is 
the maximum 3D density. Even for n 3 pj ~ 10 15 cm~ 3 , 
the lifetime r can greatly exceed seconds when approach- 
ing the TG regime. For example, this is the case for 87 Rb 
(a 3 £> ~ 10 -29 cm 6 /s) optically trapped with uj w 100 
kHz (l « 200 A), assuming L w 100 fim and N = 200. 
Then one has 7 w 10 and Eq. I jlip predicts a reduction 
of the three-body rate by more than three orders of mag- 
nitude. 

We now turn to phase coherence of a ID Bose gas in 
a harmonic potential V(x) — mu> 2 x 2 /2. We consider 
the case of T = and rely on the hydrodynamic ap- 
proach [lij l in which long- wave properties of the ID fluid 
are described in terms of two conjugate fields, density 
fluctuations Sn and phase <fi. They satisfy the commuta- 
tion relation [8n(x], exp (i(f>(x'))] = 8(x — x') exp (i(fi(x)). 
We assume the Thomas-Fermi regime and use the local 
density approximation [l?l l2f)j|: the mean density n(x) 
is obtained from the local equation of state fj,(n(x)) = 
A*o — V(x), where fx(n) is the chemical potential for 
the Lieb-Liniger problem. The density is non-zero only 
within the Thomas- Fermi radius Rt f = \j2nQ/m(jj 2 , and 
the normalization condition n(x)dx = N gives a 

relation between hq and N. Equations of motion for the 
fields 5n and <p follow from the quantum Hamiltonian: 



H i=2^. / dx K(*) 2 + vj{d x 4>) 2 ) = /'2_, <>.//,//.. 
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where Qj and bj are frequencies and annihilation opera- 
tors of elementary excitations. The quantities vn(x) — 
(irh)^ 1 d /i / ' dn and vj(x) = Trhn(x)/m determine the lo- 
cal sound velocity vm(x)vj{x) and the local Luttinger 
parameter K{x) — y/ v,j{x)/vn(x). The Hamiltonian H q 
is a generalization of the effective harmonic Hamiltonian 
of Ref. |l4| to a non-uniform system. 

Using the density-phase representation for the field 
operators, we calculate the one-particle density matrix 
gx(x,x') = (* t (x)*(x')) for \x-x'\ > l c . As the density 
fluctuations are small, this matrix reduces to 



9l (x,x') = y/n(x)n(x')exp{-((<f>{x) - 4>{x')f) /2}. 

The phase operator is given by its expansion in eigen- 
modes labeled by an integer quantum number j > 0: 



(12) 



where we have introduced a dimensionless coordinate y = 
x/Rtf- The eigenfunctions fj are normalized by the 



4 



1 



1 ^"""""""^^ 

\ 


= 0.1 




= 1 ^\ 


i i i i 


= 10 ^^^1 



0.2 0.4 0.6 0.8 



y = x/ R T f 

FIG. 2: The density matrix gi(x, —x) for N = 10 4 and various 
values of 70. The solid curves show numerical results, and the 
dashed curves the results of the quasiclassical approach. 

condition $_ % dy (v N (0)/v N (y)) fj(y)fj'(y) = %' . From 
the Hamiltonian H q we obtain the continuity and Euler 
equations which lead to the eigenmode equation, 

(1 - y 2 )/j' - (2y/f3(y))f> + (2f7 2 //%V 2 )/ 3 - 0. (13) 

The quantity (3{y) = dhifi/dhin is determined by the 
local parameter 7(2;) = mg/h 2 n(x). In the TG regime, 
we have (3 — 2, and in the GP regime (3=1. The 
coordinate dependence of (3 is smooth, and we simplify 
Eq. (fT~S|> by setting (3{y) = (3 , where /3 is the value 
of (3 at maximum density. This simplification has been 
used [2(], H3 to study the excitation spectrum of trapped 
ID Bose gases. Then Eq. ()lrt|) yields the spectrum 
O 2 = LU 2 (j(3o/2)(j + 2/Po — 1), and the eigenfunctions 

fj(y) are Jacobi polynomials Pj a ' a \y) with a = 1/(3q — 1. 

Using Eq. ljl"2*|) . the mean square fluctuations of the 
phase ((</>(£) — 4>(x')) 2 ) are reduced to the sum over j- 
dependent terms containing eigenfunctions fj and eigen- 
frequencies fij. For the vacuum phase fluctuations this 
sum is logarithmically divergent at large j , which is sim- 
ilar to the high-momentum divergence in the uniform 
case. Accordingly, we introduce a cut-off j rnax following 
from the condition hftj ps min{/i(a;), /it(x')} and ensuring 
a phonon-like character of excitations at distances x and 
x' . The vacuum phase fluctuations have been calculated 
by using two approaches: numerical summation over the 
eigenmodes with fj , flj from the simplified Eq. (|13[) , and 
quasiclassical approach assuming that the main contri- 
bution comes from excitations with j 1. In the lat- 
ter case, for x' = —x we obtain — </>(— x)) 2 ) w 
K^ 1 (x) h\\\2x\ 1 1 r jx)\ : which is close to Haldane's result 
for a uniform system |l4| with the Luttinger parameter 
K{x) and correlation length l c {x). 

The dependence of g\ on the dimcnsionlcss coordinate 
y is governed by two parameters: 70 = 7(0) and the 
number of particles N. In Fig. we present the quan- 
tity gi(y, —y)/n(y) for N = 10 4 and various values of 70. 
As expected, the phase coherence is completely lost in 
the TG regime (70 ^> 1). Moreover, on a distance scale 



x r~j Rtf the coherence is already lost for 70 ~ 1. Ther- 
mal fluctuations of the phase are readily included in our 
scheme and will be discussed elsewhere. 

In conclusion, we have found an enhanced stability of 
a trapped ID Bose gas in the Tonks- Girardeau and in- 
termediate regimes and described the reduction of phase 
coherence in these regimes. 
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